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FIELD OF THE INVENTION 

The present invention relates to a method of evaluating a capillary pressvire curve of 
the rocks of an underground reservoir from measurements on rock debris taken 
therefrom. 

BAKGROUND OF THE INVENTION 
Laboratory measurements on cores or cuttings 

Measurement of petrophysical parameters such as the permeability, the porosity and 
the capillary properties on rock fragments obtained while drilling a well through an 
imderground formation is an interesting opportunity for operator companies to rapidly 
obtain a first petrophysical characterization of producing zones traversed by the well. 

Patent FR-2,809,821 filed by the applicant describes a system of evaluating 
physical parameters such as the absolute permeabiUty of porous rocks of an 
underground reservoir zone, from cuttings retumed to the surface in the drilling mud. In 
an enclosure where the cuttings are dipped in a viscous fluid, some of this fluid is 
injected at a pressure that increases with time until a predetermined pressure threshold is 
reached, so as to compress the gas trapped in the pores of the rock. This injection stage 
is followed by a relaxation stage where injection is stopped. The pressure evolution 
during the injection process being modelled from initial values selected for the physical 
parameters of the cuttings, a computer adjusts them iteratively so as to obtain the best 
possible agreement between the modelled pressure curve and the pressure curve really 
measured. 
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Patent application FR-02/0023 filed by the applicant describes another method of 
^ - . evaluating physical parameters such as the absolute permeability and the porosity of the 

rocks of an underground reservoir zone, also from cuttings. An enclosure containing the 
rock fragments and filled with a viscous fluid is communicated with a vessel containiag 
5 the same fluid at a predetermined pressure so as to compress the gas trapped in the 
pores of the rock. The time of application of this pressure, according to whether it is 
short or long, allows to measure either the pressure variation in the enclosure or the 
variation of the volume actually absorbed by the rock fragments. The pressure or 
volume evolution in the enclosure is then modelled from initial values selected for the 
10 physical parameters of the fragments, and the values of the physical parameters of the 
rock fragments are iteratively adjusted so that the modelled evolution best adjusts to the 
measured evolution of the physical parameter in the enclosure. 

In the sphere of petrophysical characterization, the capillary pressure is also a very 
important datum fopoperators because it conditions : 

15 - the initial distribution of the fluids in the reservoir from the aquifer zone (referred to 
as woe, for water-oil contact, by the man skilled in the art) to the upper part of the 
reservoir (transition zone). According to the capillary pressure curve associated with a 
reservoir rock and the nature of the fluids in place, this transition zone can extend over 
some meters or some ten meters, which has an important effect on the determination of 

20 the accumulations in place, 

- the input pressure of a rock, which is particularly important for cap rocks. For 
example, for a gas storage tank, the input pressure of the cap rocks conditions directly 
the allowable overpressure in the storage levels without having leaks. 
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With the current techniques, the capillary pressure curve is obtained by means of 
laboratory measurements on reservoir cores. These methods are expensive because of 
the coring operations as well as the measurements on the cores, and the results are often 
available only several months after drilling. 

5 Approaches to rapidly obtain the capillary pressure curve 

However, there are alternative methods described in the literature for evaluating the 
capillary pressure curve rapidly, either during drilling or sUghtly later. 

The most commonly used approach consists in using the mercury porosimetry 
technique for measuring the air/mercury capillary pressure curve Pc directly from 
10 cuttings. However, the curve obtained is significantly different from the reference curve 
obtained on cores with high wetting fluid saturations. Besides, this approach is based on 
the use of mercury, which is extremely polluting and progressively forbidden by the law 
in many countries, \^ch poses a major problem for applying this technique in the near 
future. 

15 Another known method uses the Nuclear Magnetic Resonance (NMR) technique to 

rapidly estimate the capillary pressure curve from stratigraphic data measured in the 
well shortly after drilling. It is notably described in the following publications : 

- Bowers, M., A. et al. : « Prediction of permeability from capillary pressure curves 
derived with NMR », 17 September 1998, 

20 - Marshall, D., et al. : « Method for correlating NMR relaxometry and merciuy 
injection data », SCANo. Society of Core Analysts International Symposium 1995, 
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- Volokitin, Y., W. J. et al. : « A practical approach to obtain 1^ drainage capillary 
pressure curves from NMR core and log data », SCA No. Society of Core Analysts 
International Symposium 1999. 

The NMR relaxation signal is first converted in terms of pore size distribution, then 
5 in terms of threshold size distribution, which allows to calculate a pseudo-capillary 
pressure curve. This approach has been tested on several samples of known curve Pc. 
The results show that a good agreement with the reference curves can only be obtained 
by means of a rigorous calibration stage to be carried out case by case according to the 
nature of the rocks studied. This calibration stage is necessary owing to the uncertainty 
10 on : 

- the NMR signal-pore size distribution conversion which depends on the value of the 
surface relaxivity which is variable according to the rocks, and 

- the pore size distribution-threshold size distribution conversion which depends on 
the nature of the rock and on the diagenesis process. 

15 This approach is therefore not recommended in a predictive exploration context. In 

any case, it would not be applicable to cuttings. 

Image analysis has also been the subject of work intended to obtain a curve Pc. The 
porous medium is first prepared in form of a thin section photographed by scanning 
electron microscopy or SEM. The image obtained is then analysed so as to determine 
20 parameters representative of the proportion and of the shape of the voids in relation to 
the rock. In particular, it is possible to determine a threshold size distribution to 
reconstruct a pseudo-capillary pressxu'e curve Pc. The main limitation of this method is 
the two-dimensional (2D) nature of the thin section, whereas the capillary pressure is by 
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definition a three-dimensional (3D) property. Besides, tihis technique requires quite 
heavy conditioning, which is not really compatible with a result obtained slightly later. 
Image analysis could be applicable on cuttings but it would require carefiil calibration to 
acquire a good predictability. 

5 Finally, it can be noted that the centrifuging technique is sometimes applied in the 

field to cuttings, but in order to extract the largest possible amount of drilling fluid firom 
the rock to minimize pollutant discharges to the environment and to limit the cost by 
recycling the drilling fluid recovered. As far as we know, no centrifugation on cuttings 
has been considered in order to detemune capillary properties. 

10 SUMMARY OF THE INVENTION 

The method according to the invention allows to determiae the capillary pressure 
curve of rocks of an imderground reservoir from measurements on rock debris or 
fragments (such as cuttings) taken therefrom, over the total saturation range of these 
rocks, within a short period and at a low cost, from these measurements. It comprises : 

15 - measuring the permeability k of the rock debris, 

- measuring the capillary pressure curve Pc as a ftinction of the saturation of the rock 
debris iaitially saturated with a fluid by subjectiag them to centrifugation, and 

- parametrizing a capillary pressure curve Pc satisfying empirical relations depending 
on adjustable parameters, that is constrained to adjust to an asymptotic part of the 

20 capillary curve measured by centrifugation, and to the value of permeability k measured 
on the cuttings, so as to obtain the whole of the capillary pressure curve. 



ParametFizing the curve is advantageously carried out by selecting by default a set 
of said parameters allowing calibration on the asymptotic part of the capillary pressure 
Pc with low saturations, and by modifying the parameters step by step so that the 
estimation of the permeability given by one of the empirical relations used is best 
adjusted with the measurements of permeability k carried out on rock debris and with 
this asymptotic part. 

Permeability k of the cuttings is measured for example from measurements of the 
pressure variations in a vessel filled with a fluid containing the cuttings after it has been 
communicated for a predetermined period of time with a tank containing the same fluid 
under pressure, and from the volume actually absorbed by the cuttings, and from 
modelling the evolution of the pressure or of the volume in the vessel, from initial 
values selected for the physical parameters of the cuttings, which are iteratively adjusted 
so that the modelled pressure evolution best adjusts with the measured evolution of the 
physical parameters of the cuttings. 

The method is notably advantageous in that it provides the capillary pressure of the 
rocks on the basis of simple cuttings that are more readily available and less e3q)ensive 
to obtain. The results are also obtained much more rapidly than vnth cores. 

BRIEF DESCRIPTION OF THE FIGURES 

Other features and advantages of the method according to the invention will be 
clear from reading the description hereafter of an embodimeiit given by way of non 
limitative example, with reference to the accompanying drawings wherein : 

- Figure 1 shows the evolution of the NMR signal during various centrifugation 
stages carried out on cuttings. 
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- Figure 2a shows the curves Pc ohtained on a single rock (B7) from cuttings and 
from a core, 

- Figure 2b shows the conq)arison between the reference curve Pc measured on a core 
and the reconstructed curve Pc over the total saturation range from the centrifugation 

5 measurements and the permeabihty measurement, 

- Figures 3a and 3b show results comparable to those shown in Figs.2a and 2b 
respectively for another rock GDVl, 

- Figures 4a and 4b show results comparable to those shown in Figs.2a and 2b 
respectively for another rock Rotl, 

10 - Figures 5a and 5b show results comparable to those shown in Figs.2a and 2b 
respectively for another rock St Max, 

- Figure 6 shows the comparison between the reference permeabilities and the 
permeabilities measured on cuttings within the scope of the method described in the 
aforementioned patent appUcation FR-02/0023, and 

15 - Figure 7 diagrammatically shows the flowchart for urqplementing the method. 

DETAILED DESCRIPTION 

The method intended for fest evaluation of a capillary pressure curve Pc from rock 
fragments or cuttings according to the invention is illustrated by Figure 7. The method 
is based on two experimental measuring stages, followed by a parametrizing stage by 
20 reference with known curves, to reconstruct curve Pc over the total saturation range. 
During one of the measuring stages, capillary pressure data Pc are acquired by 
centrifuging the initially water-saturated rock fragments. The other experimental stage 
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allows to calculate the permeability value k of the rock from the method described in 
the aforementioned patent application FR-02/0023. Reconstruction of curve Pc over the 
total saturation range is carried out by using a parametrized form. The parameters of the 
curve are determined in such a way that the curve matches the capillary pressiire data Pc 
obtained e3q)erimentally and the estiroated permeability value from the curve Pc 
obtained according to a known method such as Thomeer's or Swanson-Kamath's 
method, which are reminded below, with the value measured on the rock fragments. 
Several appUcation cases are presented, with show the very good agreement obtained 
with reference curves without any particular previous calibration procedure. 

T) Measurement of Pc by centrifugation from the cuttings 

To implement the method, it is possible to use standard centrifugation means or 
more sophisticated means with automatic monitoring of the volumes of fluid produced, 
such as those described for example in patents EP-603,040 (US-5,463,894), FR- 
2,763,690, FR-2,772,477 (TJS-6, 185,985) or FR-2,798,734 filed by the applicant. 

The cuttings \vlnch have come up to the sur&ce during the drilling operation are 
first cleaned with solvents in a Soxhlet type device, then dried and saturated with 30 g/1 
brine. The cuttings are then drained in a damp cloth so as to remove the water trapped 
between the various cuttings d, then placed in a cell or cup fastened to the end of a 
rotating arm. The water expelled by centrifugation from the cuttings flows throxigh a 
grate and it is collected at the base of the cup. The experimental data are acquired in the 
same way as in the context of centrifugation on a core. For a centrifiigation stage (given 
rotating speed), the evolution of the water production is meastired until no significant 
variation can be observed any longer, then the rotating speed is increased to start a new 
stage. 
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As can be seen in Figure 1, a progressive decrease of the NMR signal and a shift to 
the short relaxation times T2 are observed, which express a desaturation of the porous 
medium with the increase of the rotating speed. This type of measurement shows that 
there definitely is a capillary contact between the cuttings, which allows to measure Pc 
5 by centrifugation from cuttings. 

The volume of water produced during the ejqjeriment is converted to saturation data 
from the volume of water initially contained in the cuttings. The latter is determined by 
weighing (diflference in the weight of the cuttings before and after saturation) or directly 
by NMR measurement. 

10 Figures 2a to 5a show the result of experiments carried out from model cuttings, 1 

to 2 mm in size, manufactured in the laboratory from rocks of known properties for 
which a curve Pc conventionally measured by centrifugation of a core is available. It 
can be seen that a good agreement is obtained with the reference curve at the level of 
the asymptotic part (low wetting flxiid saturation). On the other hand, a big difference is 

15 observed for higher wetting fluid saturations. Results equivalent to those obtained 
within the context of porosimetry measurements using mercury on cuttings are thus 
obtained, without pollution risks. 

The measured capillary pressure curve Pc is however representative only on the 
asymptotic part. A reconstruction procedure is therefore necessary to evaluate the 
20 behaviour of the curve over the total saturation range. 
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II) Measurement of the permeability k of the cuttings 

The method described in the aforementioned patent application 02/02,242 is applied 
to measure the permeability of the cuttings. The cuttings are therefore dipped in a 
containment enclosure containing a viscous fluid. The enclosure is then communicated 
5 with a vessel containing the same fluid under pressure, so as to compress the gas 
trapped in the pores of the rock. According to a first embodiment, this communication 
period can be very short and followed, after a latency time, by the measiirement of the 
pressure evolution in the enclosure. According to another embodiment, the 
commimication period can be long enough to allow to observe and measure the 
10 variation of the volume actually absorbed by the cuttings. 

The evolution of the pressure or of the volume in the enclosure is then modelled 
from initial values selected for the physical parameters of the cuttings, and the values of 
the physical parameters of the cuttings are iteratively adjusted so that the modelled 
evolution best adjusts with the measured evolution of the physical parameter in the 
15 enclosure. 

This procedure gives excellent results. The permeability values k of the cuttings are 
totally in accordance with the reference measurements obtained from cores. 

III) Reconstruction of the total capillary pressure curve 

In this third stage, we synthesize the previous measurements that have allowed to 
20 construct the asymptotic part of the capillary pressure curve Pc and the measurements of 
permeability k using empirical relations reputed to model the physical parameters of 
rocks well. The following publications : 



«' 
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Thomeer, J.H.M. : 'Introduction of a pore geometrical factor defined by the 



- Thomeer, J.H.M. : "Air permeability as a function of three pore network 
parameters". Trans AIME, vol. April, pp.809-814, 1983 

5 describe methods for evaluating the permeability firom a capillary pressure curve. 
Capillary pressure curve Pc is modelled m the following form : 



where : 

G is a shape parameter for taking accoimt of the curvature of the capillary pressure 
10 curve (related to the shape of the pore size distribution), 

- Pd the displacement pressure extrapolated to Sug equal to zero, and 

- Vboo the percentage of volimae occupied by the mercury at the end of the experiment 
at an infinite capillary pressure (equal to <ti x S^g), 

The three parameters of the model being related to the permeability by the 
15 following expression : 



capillary pressure curve''. Trans AIME, vol. March, pp.73-77, 1960, and 



i>,=P^xexp(- 



ifc = 3.8068 xG 



1-1.334 



X 




In the following publication : 

- Swanson, B.F. : "A simple correlation between permeability and mercury capillary 
pressures", JPT, vol. December, pp.2498-2504, 1981, 



12 



the author proposes correlating the value of the permeability with the maximum value 
of ratio (Vi/Pc) on the mercury porosimetry curve. This particular point generally 
corresponds to the regime change that occurs at the end of the percolation regime, just 
before the significant capillary pressure increase. As regards the capacity of fluid flow 
in the porous medium, this point is particularly important because it represents the pore 
size for which the entire pore network is connected and which therefore controls the 
flow. The most general correlation provided by the author is given by the expression 
(Swanson 1981) : 



The following publication : 

- Kamath, J. : '^Evaluation of accuracy of estimating air permeability from mercury 
injection data", SPE Formation evaluation, vol.7,4, pp.304-310, 1992 

also relates to the comparative evaluation of the approaches using empirical correlation 
or physical models to determine the permeability value from a mercury porosimetry 
curve. The best eigreement is obtained with a new correlation based on the 
aforementioned characteristic length by Swanson 1981, mentioned above : 




* = 413xZif 



ifk< ImD 
if k> ImD^ 



Lnax being defined as foUows : 




and 




100-5, 
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where : 

- k : exponent expressing the curvature of the capillary pressure curve (related to the 
shape of the pore size distribution), 

- : the displacement pressure extrapolated to Si^ equal to zero, and 
5 - Sr : residual saturation occupied by the wetting fluid (%). 

To parametrize the empirical capillary pressure curve Pc obtained from Thomeer's, 
Swanson's or Kamath's approaches, it is constrained to adjust to the asymptotic part 
obtained by centrifiigation during the first stage of the method. The entire curve is 
constrained using also the value of permeability k measured on cuttings duriDg the 
10 second stage of the method, which is compared with the result of the empirical 
relations. The parameters of the capillary pressure Pc are then modified until both the 
measured asymptotic behaviour and the permeability estin^tion are met, which allows 
to constrain the capillary pressure curve over the total saturation range Sw. 

The inversion process starts with a set of default parameters which allow to 
15 calibrate the asymptotic behaviour of the capillary pressure Pc with low water 
saturations. These parameters are then modified step by step (mainly the input pressure 
Pe or Pd and the shape factor X or G) so that the estimation of the permeability given by 
one of the previous relations is in good agreement with the permeability measurement 
obtained on cuttings while keeping a good agreement with the measurements of Pc at 
20 low water saturations. 
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1.1 Results obtained 

Figures 2b, 3b, 4b and 5b show the comparison between the curve Pc reconstructed 
according to the previous procedure and the reference curve Pc obtained on a core. It 
can be seen that, whatever the example considered, the reconstruction method allows to 
obtain a pertinent evolution of Pc over the total saturation range and in particular at high 
water saturations, whatever the permeability of the rock. 
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